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Ordered phases such as charge- and spin-density wave state accompany either full or partial gapping of Fermi 
surface (FS) leading a metal-insulator or metal-metal transition (MMT). However, there are examples of MMT 
without any signatures of symmetry breaking. One example is Na2Ti 2 St)20, where a partial gapping of FS is 
observed but a density wave ordering has not been found. Here we propose a microscopic mechanism of such 
a MMT which occurs due to a momentum dependent spin-orbit coupled molecular orbital polarization. Since a 
molecular d orbital polarization is present due to a small spin-orbit coupling of Ti, there is no spontaneous sym¬ 
metry breaking involved. However, a sharp increase of polarization happens above a critical electron interaction 
which gaps out the d orbtial FS and reduces the density of states significantly, while the rest of FS associated 
with Sb p orbtials is almost intact across MMT. Experimental implications to test our proposal and applications 
to other systems are also discussed. 

PACS numbers: 71.20.Be, 71.70.Ej, 75.30.Gw, 75.70.Tj 


Introduction - Fermi surface (FS) is one of most impor¬ 
tant concepts in solid state physics, and exists in every metal, 
semimetal, and doped semiconductors. Many ordered phases 
such as charge- and spin-density waves, and superconductiv¬ 
ity can be regarded as FS instabilities. Such ordered phases 
accompany broken symmetries involving either translation, 
time-reversal, or charge conservation symmetries. As a conse¬ 
quence, either full or partial gapping of FS occurs with anoma¬ 
lies in spin susceptibility, speicific heats, and resistivity at a 
critial temperature T c [ ]. 

However, there are examples of metallic systems which un¬ 
dergo a phase transition to a metallic state with a partial gap¬ 
ping of FS at lower temperatures without any detection of 
spontaneous symmetry breaking. A widely studied example is 
URu 2 Si 2 where both / and d orbitals are relevant and a low- 
temperature metallic phase is associated with a hidden order, 
implying difficulties of identifying the order parameter[2, 3]. 
A less studied material is Na 2 Ti 2 Sb 2 0 (NTSO), where Ti 
d orbitals and Sb p orbtials play a major role in determin¬ 
ing physical properties[4, 5]. NTSO shows a metal-metal 
transition (MMT) around T c ~ 115 K, where the temper¬ 
ature dependence of the resistivity above and below T c in¬ 
dicates its metallic behavior, while the increase of resistiv¬ 
ity and reduction of susceptibility at T c imply the reduction 
of density of states (DOS) at the Fermi level[6-8]. Specific 
heat versus T shows a sharp peak at T c , signifying substan¬ 
tial entropy change[ ( ]. Previous studies have proposed pos¬ 
sible charge/spin-density wave instabilities for the origin of 
MMT [10-12], however there is no experimental evidence of 
charge/spin-density wave order. 23 Na nuclear magnetic reso¬ 
nance measurement revealed no sign of enhanced spin flcuta- 
tions or magnetic order[13], and neutron and x-ray diffraction 
data found only changes in lattice constants across T c [ 6, 14]. 

In this work, we provide a microscopic route to MMT with¬ 
out spontaneous symmetry breaking in NTSO [15]. Fig. 1 
illustrates our proposal; a momentum-dependent polarization 
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FIG. 1. (Color online) Schematic figures illustrating evolution of 
electronic structures in NTSO upon inclusion of spin-orbit cou¬ 
pling (SOC) and on-site Coulomb interaction U. (a) and (b) show 
schematic band structures near three special k points and Fermi sur¬ 
faces without and with the presence of SOC. As depicted in (c), the 
Coulomb interaction enhances the polarization ((p)) of the spin-orbit 
coupled molecular orbital (SO-MO) states near M point. Line colors 
depict the orbital character as shown in the figure, and thicker lines 
represent the states affected by inclusion of SOC. 


of molecular orbitals (MO), which is induced by the coop¬ 
eration of spin-orbit coupling (SOC) and the on-site Coulomb 
interaction inherent in Ti 3d orbital, gaps out a significant frac¬ 
tion of FS while leaving the rest of FS intact, so that the sys¬ 
tem remains metallic. As shown in Fig. 1(a), the system has 
three Fermi surfaces with atomic Sb p, Ti d , and MO originat¬ 
ing from Ti d orbitals at T, X, and M points in the Brillouin 
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zone, respectively. Degeneracies at T and M points are pro¬ 
tected by crystal and SU(2) spin rotational symmetries, thus 
introduction of SOC of Sb p and Ti d lifts degeneracies at T 
and M respectively. Due to the smallness of Ti SOC, it makes 
a tiny gap at M point and the MO states evolve into spin-orbit 
coupled MO (SO-MO) states as shown in Fig. 1(b). Inclu¬ 
sion of the Coulomb interaction sharply enhances the SO-MO 
polarization, such that a significant gap opens at M point, re¬ 
ducing the FS area in T and X points to balance the charge 
as shown in Fig. 1(c). Neither lattice nor time-reversal sym¬ 
metry (TRS) are broken in the resulting metallic state. Due to 
the lack of spontaneously broken symmetry MMT is a first- 
order phase transition, which we will discuss later. The spin- 
orbit coupled and anisotropic nature of the SO-MO states can 
explain the strong anisotropies in electronic and magnetic re¬ 
sponses observed in recent experiments [9, 14]. While our the¬ 
ory is applied to NTSO, it can be generalized to a system with 
different orbital characters which compose FS, and undergoes 
a partial gapping of FS across MMT. 

Structure and ab-initio calculations - Crystal structure of 
NTSO consists of neighboring Ti2Sb20 layers with Na ions 
intercalated in between. As shown in Fig. 2(a), each Ti 2 Sb20 
layer is composed of a T^O square lattice and two TiSb 2 rib¬ 
bons perpendicular to each other, with the unit cell contain¬ 
ing two Ti sites in a NTSO layer (Til and Ti2 shown in Fig. 
2). The space group for the whole unit cell, which contains 
two NTSO layer, is IA/mmm (No. 139), and for an isolated 
NTSO layer the layer group is PA/mmm (No. 123). 

Previous ab-initio studies showed that, strong hybridization 
between the Ti d and Sb p orbitals yields dispersive bands 
for both states so that the system has multiple Fermi surfaces 
with d and p orbital character as shown in Fig. 1 [10, 1]]. 
Contrary to the p orbital pocket at T, which shows three- 
dimensional shape, the d orbital pockets show weak disper¬ 
sion along the layer-normal direction and considered as quasi- 
two-dimensional Fermi surfaces. While we simplify the sys¬ 
tem by choosing an isolated NTSO layer as a unit cell, com¬ 
parison between the band structures from the full and our two- 
dimensional unit cells in Supplementary Material [16] shows 
almost unaffected d bands near the zone boundary by the layer 
stacking. 

Fig. 2(a) shows the Ti d orbitals which contribute 
to the bands near the Fermi level; {\d xy )i, \d Cg )i} and 
{\d xy ) 2 : |^e g ) 2 } for Til and Ti2 respectively, where the sub¬ 
scripts 1 and 2 denote the Ti atoms to which the orbitals be¬ 
long. The hybrid \d eg )i^ orbitals are defined as linear com¬ 
binations of e g orbitals such that \d e )i = —a\d x 2 _ y 2 )i — 
b\d 3 z 2 _ r 2) 1 and \d eg ) 2 = +a\d x 2 _ y 2) 2 - b\d 3 z 2 _ r 2 ) 2 , where 
the real coefficients a and b are determined by the ligand 
fields. Note that, |cZ e ^) 1,2 are dominated by \d 3 y 2 _ r 2) 1 and 
\d 3 x 2 _ r 2) 2 respectively in our system, and \d xy )i and \d e )i 
are transformed to — \d xy ) 2 and \d eg ) 2 respectively by the C 4 
rotation depicted in the figure. Left panel of Fig. 2(b) shows 
the bands without SOC where the orbital weights of Sb p and 
\de g )i ,2 states are represented as the size of the corresponding 
symbols. The \d xy )i^ 2 weight is concentrated on the X-M line 



FIG. 2. (Color online) (a) Crystal structure of single Ti2Sb20 layer 
with dominant Ti d orbitals near the Fermi level — {| d xy ) 1 , | d eg ) 1 } 
and {\d X y} 2 , \d eg } 2 } for Til and 2, respectively — depicted in the 
figure. Note that Til and Ti2 sites are transformed to each other by 
the C 4 rotation along the z axis, (b) Band structure of single layer 
NTSO without including SOC and Coulomb interaction. Left panel 
shows orbital weights of Sb p and Ti \d eg ) 1,2 within the Bloch states 
with the weight proportional to the size of corresponding symbols. 
Right panel shows Ti | d xy )\p orbital weights, (c) Schematic figures 
of two degenerate a and states at M point depicted in the y/2 x \[2 
times enlarged unit cell (original and enlarged unit cells represented 
as dashed and dotted square, respectively) as MOs. Bands with the 
weights of the a and f5 MO states are shown in (d). 


as shown in the right panel in the figure. 

Both the T and M points have D 4 ^ point group symmetry 
which allows the presence of quadratic band touching points 
located just below the Fermi level, as can be seen in Fig. 2(b). 
At T point, the band touching consists of Sb atomic p x and p y 
orbitals where degeneracy of the orbitals is compatible with 
the Sb site symmtery (C± v ). On the contrary, the band touch¬ 
ing at M point has dominant d character and described by 
molecular orbitals (MO). Without considering spin, two de¬ 
generate Bloch states at M point which we denote as | a) and 
\/3) are expressed as follows, 

\a) « y> kM ' R (\d xy )n,i + |d efl )R, 2 ) , 

R 

1/3) PS y^ e *k M -R ( \de g )r,1 - \dxy) R., 2 ) , 

R 

where 1 cm = (^, ^), R and { 1 , 2 } are indices for Bravais 
lattice and Ti sublattices respectively. Fig. 2(c) shows the 
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FIG. 3. (Color online) (a) Band structures and PDOS projected onto 
the SO-MO states without (left panel) and with including the on-site 
Coulomb interaction (U e ff = 3 eV, right panel) in the presence of 
SOC. Orbital weights of Sb j = 3/2, 1/2, and SO-MO states are 
depicted as the size of the symbols indicating each state. Magnified 
view of bands near M point is shown in the inset, (b) A plot of nor¬ 
malized MO polarization P between the SO-MO states as a function 
of U e ff • (c) Fermi surfaces of systems without (left panel) and with 
the Coulomb interactions (right panel). 


schematic illustrations of the states, where the dashed and 
dotted squares depict the primitive and the enlarged unit cell. 
Like the Sb p x and p y orbitals at T point, they belong to the E u 
irreducible representation, which are odd under spatial inver¬ 
sion and transforms as (C^) 2 \a) = C 4 1/3> = — | a). There is 
additional degeneracy due to the SU(2) symmetry in the spin 
subspace without the presence of SOC, so the band touching 
at M point is fourfold degenerate and protected by the SU(2) 
and crystal symmetries. Unless either the spin or the crystal 
symmetry is lifted, the degeneracy at M point remain robust, 
and our LDA+U calculation without including SOC and keep¬ 
ing paramagnetic constraint confirms it (See Supplementary 
Material). 

Degeneracy lifting at T and M points via SOC - Below we 
discuss the effect of SOC. Intuitively, one expects a relatively 
large impact of SOC on Sb p orbitals while negligible on Ti 
d orbtials, as SOC in Ti is only about 20 meV. Indeed degen¬ 
erate bands near T point made of p orbitals are split, and the 
gap between the two is about 0.4 eV. Each bands are char¬ 
acterized by total angular momenta j = 1/2 and 3/2 made 
of | l z = ± 1 ) = | p x ± ip y ) and spin- 1/2 due to ligand field 
spitting of layer structure. As SOC is introduced, the SU(2) 
symmetry in the spin space is lifted, and the fourfold degen¬ 
eracy at T is split into j = 1/2 and 3/2 doublets as shown 
in Fig. 3(a). Hereafter we denote the orbital and total angu¬ 
lar momenta for the atomic and MO states as lowercase and 


capital letters respectively. 

Similarly, at M point, quenching of planar orbital moment 
components happens in the MO space {\a),\f3)}. Projecting 
the Ti total angular momentum operators L = t 111 ® fl 11 ' ® 
F' 2 ® t Tl2 ' onto the space yields L z as the only nonvanish¬ 
ing component, which is diagonalized with the basis choice 
of I Lf) s I a =p i/3) where Lf ~ ±\/3 in this system [17]. 
In the presence of Ti SOC, the | Lf) states are split into two 
spin-orbit coupled MO (SO-MO) doublets which are charac¬ 
terized by total angular momenta J — (\Lf | + |) = J + and 
(\Lf\ — |) = J~. Note that both J + and J~ form doublets 
with J+ = ± (\Lf\ + f) and J~ = ± (| Lf\ - |), respec- 
tively, and J + doublet is higher in energy than J~ by 20 meV 
at M point as shown in the inset of Fig. 3(a). 

Molecular orbital polarization enhanced by on-site 
Coulomb interaction - The size of the gap opening can be 
quantified by defining the MO polarization operator p = 
fij- — hj +, where hj± are number operators for the J ± sub¬ 
spaces. SOC in the MO space can be rewritten by employing 
the MO polarization operator as Hso = — Axip/2. Due to the 
minus sign, H$o favors positive MO polarization, but its mag¬ 
nitude is tiny because of the small Axi ~ 20 meV as can be 
seen in the PDOS of Fig. 3(a). However, once the degeneracy 
at M is lifted, the size of the splitting can be further enhanced, 
and inclusion of on-site Coulomb interaction does the role in 
this case. 

Like iTso> the LDA+U correction to the Ti on-site potential 
introduced by the Coulomb interaction can be projected onto 
the MO space. It can be rewritten in terms of p as follows; [18] 


Vu 




where U e x = U — J is the effective Coulomb interaction pa¬ 
rameter, ft = hj+ + fij- is the total number operator, and 
expectation values are obtained by integrating over the Bril- 
louin zone (detailed derivation is in Supplementary material). 
Combining it with Hso, apart from the trivial constant term, 
yields 


H$o + Vu 


A T i + Ueff (p) ^ 

-o- V, 


so that the MO polarization (p) initiated by SOC can be fur¬ 
ther increased under the presence of C7 e ff. Right panel in Fig. 
3(a) shows the band structure and PDOS projected onto the 
SO-MO states with the presence of U e s = 3 eV. The gap at M 
point is allowed due to the loss of SU(2) symmetry via SOC, 
and is greatly enhanced by the inclusion of the Coulomb in¬ 
teraction. A normalized MO polarization P = (p) / (ft) as a 
function of U e & is plotted in Fig. 3(b). One can see an abrupt 
change in P near U = 2.2 eV, where the bottom of the 
(| Lf | + \) band crosses the Fermi level and the electron-like 
Fermi pocket near M point disappears. Fig. 3(c) compares the 
Fermi surfaces with and without including U e s = 3 eV. To 
compensate the removal of the M point electron pocket, the 
size of hole-like pocket near X point also decreases, resulting 
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the half-reduction of the Fermi surface area. Such reduction 
of Fermi surface area, or equivalently the reduction of DOS at 
the Fermi level, is observed in several magnetic susceptibility 
and resistivity measurements [6, 7, 9]. 

We also comment on the X pocket, which is less affected 
by the inclusion of f7 e ff than the M pocket despite its domi¬ 
nant \d eg )i^ character. Since the X pocket coexists with the 
p-originated Y pocket, the size of X pocket is determined by 
the on-site energy of ) 1,2 orbitals relative to that of Sb p. 
The V[j term, however, behaves as an effective SOC for the 
SO-MO states, and its contribution to the |d eg )i ,2 on-site en¬ 
ergy is not significant so that the X pocket remains even after 
inclusion of the Coulomb interaction. This is confirmed by 
examining the change of the on-site energies from the Wan- 
nier orbital calculations for the Ti d states. Note also that, the 
Sb j = 3/2 and 1/2 states near Y are well polarized due to the 
large Asb, so they are less affected by the inclusion of U e ff. 

Discussion and Conclusion - The notion of orbital polar¬ 
ization in transition metal compounds, which is induced by 
the Coulomb interaction in degenerate d orbitals, usually ac¬ 
companies symmetry lowering by the orbital-lattice coupling 
such as Jahn-Teller effect[ 19-23]. In metallic systems, an or¬ 
bital polarization also occurs [24], but does not accompany a 
MMT for a fixed charge filling. The SO-MO polarization in 
our work is distinguished from previous studies since it ac¬ 
companies a phase transition without spontaneous symmetry 
reduction. 

For experimental validation of our proposal, one can take 
advantage of the anisotropy in the SO-MO states. Since 
the total angular momentum of the SO-MO states are fixed 
to be perpendicular to the Ti 2 0 plane, i.e. the absence of 
the in-plane angular momentum component, one direct con- 
squence of the SO-MO formation is an anisotropic response 
to the external magnetic fields. A recent report on the mag¬ 
netic susceptibility of NTSO shows such anisotropic behav¬ 
ior, where the H/ jz susceptibility data shows drastic en¬ 
hancement compared to the H//xy result below the transi¬ 
tion temperature [9]. Such behavior is hard to understand in 
the conventional charge-density wave picture, but is consistent 
with our SO-MO polarization scenario; due to the absence of 
the L x , y components in the SO-MO space, the states should 
be more susceptible to H//z compared to the in-plane fields. 

Measuring the branching ratio in Ti L 2 and L 3 edge x-ray 
absorption spectroscopy (XAS) below T c can be another way 
to experimentally validate the MO polarization scenario [25]. 
Note that, the unoccupied SO-MO state has total angular mo¬ 
mentum of j + = (\lf\ + |), which is close to 5/2, for each 
Ti site. As a result, L 3 edge XAS is expected to show higher 
intensity more than twice that of L 2 edge, so that the ratio be¬ 
tween the L 3 and L 2 channels should be larger then the statis¬ 
tical ratio 2:1. Also, phase-sensitive tools such as resonant x- 
ray scattering may probe the presence of SO-MO states, which 
revealed the presence of the spin-orbit-entangled j e ff = 1/2 
states in several iridate compounds with strong SOC[26, 27]. 

Finally, we would like to mention the nature of the phase 
transition at T c . Specific heat versus temperature data indi¬ 


cate significant loss of entropy at the transition. Previously 
this was considered as a signature of continuous phase transi¬ 
tion by the density wave formation, but this feature can be ex¬ 
plained within our SO-MO polarization scenario which results 
in entropy change by the reduction of the Fermi surface area. 
Since this process does not spontaneously break any symme¬ 
tries in NTSO, it should be a first-order transition. We found 
a metastable phase with P ~ 0 in the regime 2.3 < C/ e ff < 4 
eV and its presence is robust independent of parameters and 
code choices (see Supplementary Material for details on the 
metastable state). Remarkably, hystresis behavior in resistiv¬ 
ity with unknown origin was reported previously[7], indicat¬ 
ing a possible first-order phase transition. 

In summary, we propose an alternative scenario for the 
MMT in NTSO based on the SO-MO polarization induced 
by SOC and the Coulomb interaction. Our picture can be 
generalized to other systems with FS with mixed orbital char¬ 
acters and spontaneous phase transition with partial gapping 
out of FS without a symmetry breaking. Also, our picture 
calls attention to the role of SOC in 3d transition metal com¬ 
pounds, which was considered insignificant in understanding 
the physics of such systems. Indeed, there are several reports 
about the role of the ostensibly small SOC, which cooperate 
with the Coulomb interaction and becomes a crucial element 
for the underlying physics[28, 29]. Further studies on the role 
of SOC in previously known transition metal systems with 
high crystal symmetries can be an interesting future subject 
in this regard. 
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